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► We provide the first-ever measurements of the Congo wetlands water volume change. ► Wetland water is dominated by local upland runoff and much less 
from mainstem. ► Differences between the Congo wetland and the Amazon floodplain are highlighted. 
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The Congo Basin is the world's third largest in size (~3.7 million km 2 ), and second only to the Amazon River 26 
in discharge (-40,200 m 3 s -1 annual average). However, the hydrological dynamics of seasonally flooded 27 
wetlands and floodplains remains poorly quantified. Here, we separate the Congo wetland into four 3°x3° 28 
regions, and use remote sensing measurements (i.e., GRACE, satellite radar altimeter, GPCP, JERS-1, SRTM, 29 
and MODIS) to estimate the amounts of water filling and draining from the Congo wetland, and to determine 30 
the source of the water. We find that the amount of water annually filling and draining the Congo wetlands is 31 
111 km 3 , which is about one-third the size of the water volumes found on the mainstem Amazon floodplain. 32 
Based on amplitude comparisons among the water volume changes and timing comparisons among their 33 
fluxes, we conclude that the local upland runoff is the main source of the Congo wetland water, not the fluvial 34 
process of river-floodplain water exchange as in the Amazon. Our hydraulic analysis using altimeter mea- 35 
surements also supports our conclusion by demonstrating that water surface elevations in the wetlands 36 
are consistently higher than the adjacent river water levels. Our research highlights differences in the hydrol- 37 
ogy and hydrodynamics between the Congo wetland and the mainstem Amazon floodplain. 38 

© 201 1 Published by Elsevier Inc. 39 
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1. Introduction 

The Congo Basin is the world’s third largest in size (-3.7 million km 2 ), 
and second only to the Amazon River in discharge (-40,200 m 3 s _ 1 an- 
nual average). The impact and connections of this hydrologic flux with 
the region’s climate, biogeochemical cycling, and terrestrial water stor- 
age, especially in wetlands, is of great importance. For example, the ex- 
tent of the differences in chemistry, seasonality, rate and volume of 
water input to the floodplain and wetland systems from upland runoff, 
direct rainfall and mainstem flooding are likely to supply substantially 
different amounts of nutrients and other solutes (Melack & Engle, 
2009). However, the hydrological dynamics of seasonally flooded 
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wetlands and floodplains remains poorly quantified through ground ob- 55 
servations, satellite observations or modeling. As a consequence, esti- 56 
mates of the magnitude of other processes driven by such dynamics, 57 
such as methane emissions from flooded wetlands that form a significant 58 
contribution to global atmospheric methane, also cannot be well estimat- 59 
ed. Given the vast size and remote location of the Congo Basin, satellite- 60 
borne observations provide the only viable approach to understanding 61 
the spatial and temporal distributions of its water balances. 62 

Recently, Alsdorf et al. (2010) have estimated the amounts of water 63 
filling and draining from the mainstem Amazon floodplain using data 64 
from the Gravity Recovery and Climate Experiment (GRACE) and 65 
other satellite measurements. They showed that the majority of water 66 
on the mainstem Amazon floodplain is derived from the river with a 67 
much less amount from local upland runoff. However, there has been 68 
no attempt to estimate the Congo wetland water storage and its flux. 69 
In this study, we use satellite-borne observations to suggest a baseline 70 
measurement of these storages and fluxes by examining 1 ) the amount 71 
of water stored and drained from the Congo wetland, and 2) whether 72 
the water comes from rivers or adjacent upland areas. 73 

We use total storage change in the form of equivalent water height 74 
(EWH) change (Wahr et al., 1998) from the GRACE measurements 75 Q2 
(Tapley et al., 2004), precipitation (P) estimates from the Global 76 
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Precipitation Climatology Project (GPCP; Adler et al., 2003), evapotrans- 
piration (ET) estimates from the Hillslope River Routing (HRR) model 
(Beighley et al., 2009), water elevation changes from Environmental 
Satellite (Envisat) altimeter measurements, and hydrological maps 
from HydroSHEDS (Lehner et al., 2008). Measurements of inundated 
area are made from a combination of ( 1 ) the Japanese Earth Resources 
Satellite- 1 (JERS-1) Synthetic Aperture Radar (SAR) mosaics developed 
by the Global Rain Forest Mapping project (GRFM), (2) the Shuttle 
Radar Topography Mission (SRTM) Digital Elevation Model (DEM), and 
(3) Moderate-resolution Imaging Spectroradiometer (MODIS) mosaics 
(Jung et al., 2010a). Unfortunately, we have no available contemporane- 
ous in situ discharge or water stage measurements. We combine these 
satellite-based measurements to: (1) estimate the wetland storage 
changes in four regions along the Congo mainstem and its major tribu- 
taries, and (2) determine whether the water comes from rivers or adja- 
cent upland areas. 

The methods presented here are improved compared to the previous 
study over the Amazon Basin (Alsdorf et al., 2010) because 1) Hydro- 
SHEDS is used to estimate the upland area that contributes directly to 
the wetland instead of using a ratio between estimates of upland area 
compared to the wetland area; 2) more realistic ET estimates are used 
instead of a single number representing the whole basin; and 3) a hy- 
draulic analysis from altimeter measurements is also presented. We 
also use a longer time span (6 years compared to 2.5 years) of GRACE 
data. 


2. Methods 

2.1. Study area 

We select four 3°x3° study regions to cover the wetlands of the 
Congo River mainstem and its major tributaries (Fig. 1). Study region 
1 includes the Ubangi River (-3800 m 3 s -1 annual discharge, Laraque 
et al. (2001)), which is the largest right-bank tributary of the Congo 
mainstem. Study region 2 includes the Sangha River (-1600 m 3 s -1 
annual discharge, Laraque et al. (2001)) and represents the majority 
of the northern tributary wetlands. Study regions 3 and 4 include east- 
ern and southern tributaries, respectively. The box size is chosen based 
on the limit of the spatial resolution of GRACE which is determined 
from the maximum degree (n max = 60) of the Stokes coefficients. 



14° 16° 18° 20° 22° 24° 26° 28° 
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Fig. 1 . Locations of four 3° x3° study regions in the Congo Basin. Background shows to- 
pography from the SRTM C-band DEM. Intersections between Envisat altimeter and the 
Congo River are indicated with 


2.2. Wetland storage changes from satellite measurements 114 

Total storage changes for a given area, AS, are a summation of the 115 
storage changes in wetlands (AS W ), rivers (AS r ), groundwater (AS g ), 116 
and soil moisture (AS sm ): 117 

AS = AS W + AS r + ASg + AS sm . (1) 

118 

Measurements from GRACE provide AS in terms of anomalies with 120 
respect to a mean total storage value. We processed the Release 4 121 
(RL04) Center for Space Research (CSR) GRACE Level 2 (L2) data 122 
product (Bettadpur, 2007) from January 2003 to December 2008. To 123 
reduce the GRACE longitudinal stripes associated with correlations 124 
among even or odd degree Stokes coefficients at resonant orders 125 
(Swenson & Wahr, 2006), decorrelation based on Duan et al. (2009) 126 
was used. We also applied smoothing using a 3-degree Gaussian filter 127 
(Guo et al., 2010). EWHs are computed at l°xl° grid spacings, and 128 
spatially averaged over each study region. Finally, total storage anom- 129 
alies are obtained by multiplying the EWHs by the box area. More de- 130 
tails on the GRACE measurements are provided in Section 3.1. 131 

The channel storage anomalies are estimated by multiplying water 132 
stage anomalies, obtained from the Envisat altimeter, with open 133 
channel areas estimated from the classification of GRFM image data 134 
(Table 1, see discussion below). The Envisat Geophysical Data Records 135 
(GDRs) contain 35-day repeat, 18-Hz data (twenty-measurements- 136 
per-frame), which corresponds to a ground spacing of approximately 137 
350 m. The GDRs include range measurements from four different 138 
retracking algorithms. In this study, we use the retracked measure- 139 
ments from the ICE-1 retracker (Bamber, 1994), which generally 140 Q3 
performs well over inland water bodies (Frappart et al., 2006; Lee 141 Q4 
et al., 2010). The water stage anomalies over the intersections 142 
between the altimeter and the open water bodies are averaged for 143 
each tributary, and are then multiplied by the corresponding channel 144 
areas. 145 

We use 2.5° x 2.5° GPCP monthly merged precipitation rates P(t) 146 
(Adler et al., 2003), and create anomalies by subtracting a linear fit, P, 147 
to the integrated sum of P(t) for each study region (see Alsdorf et al., 148 
2010 for details). The slopes of the linear-fit lines represent six-year 149 
mean precipitation values, as summarized in Table 1. The GPCP data is 150 
derived partly from infrared and microwave satellite measurements, 151 
and it should be noted that, as stated in Beighley et al. (2011), there is 152 
a discrepancy between various satellite derived precipitation datasets 153 
over the Congo Basin in terms of their magnitudes, especially in 154 
equatorial regions, which correspond to study regions 2 and 3 in 155 
this study. For ET, we use model-based estimates from HRR. It is the 156 
sum of wet canopy evaporation, dry canopy transpiration and evapora- 157 
tion from saturated soil surfaces based on the potential ET using Pen- 158 
man-Monteith indirectly through the temperature-based method of 159 
estimating its data sources (see Beighley et al., 2009, 2011 for details). 160 
The ET rates over each Pfafstetter Level 4 sub-divisions are averaged 161 
for each of the four study regions (E) (Table 1 ). This Pfafstetter discreti- 162 
zation frame- work is a natural system, based on topographic subdivision 1 63 
of the land surface and the resulting topology of the hydrographic net- 164 
work (Verdin & Verdin, 1999). Each level of discretization results in 9 165 
sub-divisions (i.e., 4 tributaries and 5 local contributing areas to the 166 


Table 1 tl.l 

Hydrologic and geomorphic characteristics of each study region. 

tl.2 



Region 1 

Region 2 

Region 3 

Region 4 

tl.3 

Upland (km 2 ) 

83,605 

42,905 

55,297 

58,587 

tl.4 

Wetland (km 2 ) 

28,052 

68,596 

56,360 

52,914 

tl.5 

Channels (km 2 ) 

1058 

3990 

502 

2766 

tl.6 

Annual P (m year -1 ) 

1.44 

1.53 

1.87 

1.71 

tl.7 

Annual ET (myear -1 ) 

0.90 

1.01 

1.06 

0.92 

tl.8 

Contributing area (km 2 ) 

121,330 

151,596 

152,789 

141,728 

tl.9 
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main channel). The ultimate number of sub-areas is Number , but 
often less than that due to lack of network resolution at higher levels. 
The Congo Basin was ultimately delineated to Pfafstetter level 4 using 
a threshold area of ~8.1 km 2 which resulted in 5498 model units (i.e., 
sub-divisions) with a median model unit drainage area of 670 km 2 
and a mean hillslope length of 5.4 km. Anomalies of P — ET are estimated 
to be: 

P-ET = (P(t)-P) x^. (2) 

The P — ET anomalies are used to estimate the runoff from the local 
uplands (Section 3.2). 

In this study, we assume that groundwater changes associated with 
the shallow water table ( AS g ) are driven by P — ET. These changes are 
assumed to be negligible beneath wetland areas that do not drain, i.e., 
the water table is assumed to be consistently at the surface in wetlands 
that contain water from year to year. P — ET varies seasonally and is 
expected to account for water table variations in the upland areas of 
each 3°x3° box. Similarly, we assume that P — ET is forcing any soil 
moisture variations (AS sm ). Thus, our estimates of P and ET are used, 
below in Section 3.2, to account for AS g and AS sm . 

The Congo interfluvial wetlands cover a variety of vegetation and 
hydrogeomorphic environments. Most of the Congo classifications 
have been developed based on vegetation type and forest density 
(e.g., Hansen et al., 2008; Laporte et al., 1995), whereas few classifica- 
tion schemes have focused on flooding in the wetlands (Bwangoy 
et al., 2009). In this study, we use the hydrogeomorphic flood classi- 
fication of Jung et al. (2010a). 

3. Results and discussions 

3.1. GRACE measurements over the Congo Basin 

The Congo River is the only major river to cross the equator twice. In 
doing so, the basin lies in both the Northern and Southern Hemisphere 
such that it receives year-round rainfall from the migration of the Inter- 
Tropical Convergence Zone (ITCZ). After the north has its wet season in 
July-September, the ITCZ moves south and the remainder of the basin 
receives large amounts of rain. Fig. 2 shows the spatial variations in 
the storage changes from the CSR GRACE data after decorrelation and 
smoothing. It can be seen that the positive anomaly in September 
2006, which is present outside of the Congo Basin, becomes stronger 
as it moves southward and into the basin. Likewise, the positive anom- 
aly observed in the southeastern part of the basin in January 2007 


becomes stronger and widely spread over the southern boundary of 207 
the Congo Basin. This spatial pattern of the storage changes is different 208 
from that over the Amazon Basin, where the strongest positive or neg- 209 
ative annual water storage anomalies are observed to be centered inside 210 
the basin (e.g., Alsdorf et al., 2010; Han et al., 2005). 211 

We also examine the basin-averaged time series of EWH anomalies 212 
obtained using four different GRACE data products (from CSR, Jet Pro- 213 
pulsion Laboratory (JPL), GeoForschungsZentrum (GFZ), Institut fur 214 
Geodasie und Geoinformation (ITG)) using equivalent decorrelation 215 
and smoothing (Fig. 3 (top)). They generally agree in terms of their an- 216 
nual increases and decreases in the time series. In addition, all of them 217 
show a drying trend until 2006, and then a sharp increase at the end 218 
of 2006. However, there are important differences in their amplitudes. 219 
There are at least 1 cm EWH differences among the GRACE products; 220 
for example, the CSR and ITG solutions differ by at least 1 cm during 221 
the last two months of 2005. If we convert this 1 cm EWH difference 222 
to streamflow by multiplying it by the basin area (3.7 million km 2 ) 223 
and dividing it by the time duration, we get about 7000 m 3 s _1 . As a 224 
comparison, this approximately corresponds to the mean annual dis- 225 
charge of the Ohio River in the United States. Moreover, there are at 226 
least 5 cm EWH differences between the CSR and JPL solutions that 227 
last about 5 months in the first half of 2008. If we again convert this to 228 
discharge, we get approximately 14,000 m 3 s -1 which corresponds to 229 
more than one-third of the Congo River mean annual discharge. It also 230 
corresponds to about three-quarters of the Mississippi River discharge. 231 
This is a significant difference: note that the Congo and Mississippi River 232 
basins are similar in size. Furthermore, the four different GRACE prod- 233 
ucts do not produce the same errors year after year. For example, in 234 
the first half of 2006, the JPL solution has generally less EWH values 235 
than the CSR solution, but in the second half of 2006 when the trough 236 
occurs, the CSR solution values are less than the JPL values. This can be 237 
widely observed every year among all of the GRACE products. Overall, 238 
the discrepancy among the GRACE products has important implications 239 
for Congo hydrology. In addition to different data processing methods 240 
and models adopted at different institutes, these disagreements may 241 
also be due, in part, to the movement of ITCZ and the consequent leak- 242 
age of strong signal from outside of the basin (e.g., the strong positive 243 
anomaly in September 2006) or from inside of the basin (e.g., the strong 244 
positive anomaly in May 2007). This leakage is due to the truncated 245 
spectral degree (e.g., n max = 60) in the GRACE gravity field solutions 246 
and to post-processing smoothing. The leakage phenomenon can 247 
occur at all scales including the finest spatial resolution possible with 248 
GRACE. 249 

Recently, in the GRACE science community, there has been an ef- 250 
fort to use global simulations of water storage variations when 251 
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Fig. 2. Monthly Equivalent Water Height (EWH) anomalies from the CSR GRACE product after decorrelation and 3-degree radius Gaussian smoothing. The Congo Basin is shown 
with a red outer boundary. Red rectangles indicate our study regions. (For interpretation of the references to color in this figure legend, the reader is referred to the web version 
of this article.) 
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Fig. 3. (top) Comparison among GRACE EWH anomalies over the entire basin after 
equivalent decorrelation and smoothing, (bottom) Comparison among EWH anomalies 
over the study regions from CSR GRACE data. 


restoring the signal loss in GRACE, which is caused by smoothing. It 
has been proposed to estimate scale factors, by comparing the unfil- 
tered model time series with the filtered model simulations, to par- 
tially correct for the signal attenuation. We have examined the 
original filtered basin-averaged time series with the scaled time se- 
ries (data courtesy, S. Swenson) from CSR GRACE data. The scale, 
computed using the Community Land Model (CLM), averaged over 
the Congo Basin is 1.2, and thus the scaled time series has a slightly 
lower amplitude than the smoothed time series. This scale is a simple 
temporally constant number, intended to depend on the statistical 
characteristic of the model-simulated storage variations. This approach 
assumes that errors in the global hydrologic model are spatially and 
temporally randomly distributed and thus do not introduce a bias in 
the scaling factor. It is also not entirely clear that a model should be 
used to correct a measurement, especially in the case of the Congo 
Basin where model errors are less well-known compared to other re- 
gions such as the United States. Moreover, the scaled EWH anomalies 
cannot resolve the issue of discrepancy among the GRACE products. 
Therefore, in this study, we do not attempt to correct the leakage 
error or restore the signal loss due to the smoothing. Rather, we treat 
the differences among the GRACE products as the error of our storage 
change estimates. 

We now compare the EWH anomalies over our four study regions 
(Fig. 3, bottom) to examine whether the 3°x3° box size is appropriate 
and if the storage changes among them can be distinguished. From 
Fig. 3 (bottom), the EWH anomalies in region 4 are clearly different 
from the other three regions in terms of their timing and amplitudes. 
For example, in 2004 region 4 has a trough in August whereas it oc- 
curred in February over the other regions. Region 4 also has the smallest 
peak in December compared to the other regions. Although the timing 


among regions 1-3 appears to agree, there are differences in anomaly 282 
amplitudes. For example, there are about 13 cm EWH differences be- 283 
tween regions 1 and 2 during July-September 2006. Converting this dif- 284 
ference to river discharge, yields about 3000m 3 s _1 . As another 285 
example, about 9 cm EWH differences, lasting about three months at 286 
the end of 2005, correspond to about 2000 m 3 s _1 of streamflow and 287 
can be observed between regions 2 and 3. In general, the amplitudes 288 
and occasionally the timing of the major peaks and troughs are different 289 
among the study regions. This distinction supports our choice of the box 290 
size and the resultant wetland volume. 291 

3.2. Wetland water volume change and its flux 292 

We observe from Fig. 4 that the total storage anomalies from GRACE, 293 
P — ET anomalies, and river storage anomalies within a given study re- 294 
gion are well timed with each other. However, in terms of amplitudes, 295 
the river channel storage anomalies are significantly less than the 296 
GRACE anomalies, which suggests that storage changes in rivers ac- 297 
count for little of the total storage anomalies (note that river anomalies 298 
are multiplied by 5, 10, or 20). The P — ET anomaly amplitudes are sig- 299 
nificantly greater than those of the rivers and typically less than the 300 
total storage anomaly amplitudes, which suggest that P — ET accounts 301 
for an important fraction of the GRACE measured total volume change. 302 
Thus we concluded that hydrological processes associated with P — ET 303 
(e.g., runoff) are significant contributors to the total storage change ob- 304 
served in each 3° x 3° study region and that in-channel fluvial processes 305 
are not significant contributors. 306 

The wetland storage anomalies have two contributors, which are 307 
(1) direct precipitation on the wetlands as well as runoff supplied 308 
to the wetlands from the surrounding uplands, and (2) water ex- 309 
changed between the wetlands and the adjacent river channels. It 310 
should be noted that the groundwater contribution to the wetland 311 
water levels is considered in the upland P — ET runoff because the 312 
groundwater is controlled by the infiltration of rainfall. The volumes 313 
of runoff from the local uplands and direct rainfall on the wetlands 314 
can be estimated by multiplying P — ET with the contributing area 315 
or with the inundated area, respectively. The land areas contributing 316 
to the wetlands are computed using the following procedures. 317 
First, flow directions from HydroSHEDS are obtained to determine 318 
flow accumulation and the associated drainage network. Next, we as- 319 
sume that major rivers have width greater than 100 m. This threshold 320 
is chosen based on the resolution of the GRFM mosaic, which is used 321 
to extract channel areas and to compute the river storage changes in 322 
Section 2.2. Based on the relationship between the channel width and 323 
the upstream drainage area (w(m) = 0.438 A u (km 2 ) 0 - 592 , Beighley 324 
et al., 2011), this 100 m threshold approximately corresponds to riv- 325 
ers with drainage areas larger than 10,000 km 2 . We remove these 326 
major rivers and their contributing areas from the flow direction 327 
grid. Thus, we distinguish the contribution of large river drainage 328 
areas from the wetland drainage areas. Then, we extract the wetland 329 
pixels for each study region using the classification map (Jung et al., 330 
2010a). Finally, we delineate the area that drains to each wetland 331 
pixel for each study region (Fig. 5 and Table 1). Essentially, the frac- 332 
tion of contributing area that is outside of each 3°x3° study region 333 
is connected with streams having a drainage area smaller than 334 
10,000 km 2 and that drain directly to a wetland pixel. To further ex- 335 
amine whether the 100 m channel width is a reasonable number to 336 
distinguish between the contributing areas that flow to the wetlands 337 
and the contributing areas that flow to the major rivers, we tracked 338 
discharges for all of the rivers which have the contributing areas larg- 339 
er than 10,000 km 2 and that flow into our study regions (red lines in 340 
Fig. 5). As summarized in Table 2, we used the HRR model to estimate 341 
these discharges during the period 2003-2008 (Beighley et al., 201 1 ). 342 
The mean annual discharge for all of the major rivers from a unit con- 343 
tributing area is estimated at 0.016 m 3 s _1 km -2 . For a drainage 344 
of 10,000 km 2 , this corresponds to 160 m 3 s _1 of discharge. So, the 345 
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Region 1 
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Region 3 


-River Storage (x20) 
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2003 



Fig. 4. Time series of satellite-based measurements of Congo hydrology for each study region. Red lines represent EWH anomalies from CSR GRACE data, and black lines are P — ET 
anomalies. Blue lines show river storage anomalies, and they are multiplied by 5, 10, or 20 for visual clarity. The river storage anomalies in this plot are generated by weighted 
averages of river stage anomalies with ratios between the channel area and the box area as the weights. (For interpretation of the references to color in this figure legend, the reader 
is referred to the web version of this article.) 


346 1 00 m river width threshold (or 1 0,000 km 2 of contributing area) leads 

347 to rivers having a discharge greater than 160 m 3 s _ 1 and which do not 

348 directly flow into the wetlands. 

349 Fig. 6 shows a comparison of water volume anomalies for the wet- 

350 lands, rivers, and local upland runoff. In each of the four plots, the am- 

351 plitudes of the river storages are negligible compared to the GRACE 

352 and P — ET anomalies. We suggest that this amplitude is not sufficient 



Fig. 5. Colored boundaries represent the contributing land area draining to the wet- 
lands in each study region with the corresponding major rivers (i.e., widths greater 
than 100 m or contributing areas larger than 10,000 km 2 ; shown as red lines) and 
their drainage areas excluded from the study regions; see Table 2 for the hydrologic 
characteristics of the 19 excluded drainages. (For interpretation of the references to 
color in this figure legend, the reader is referred to the web version of this article.) 


to supply any significant water volumes that would sufficiently account 353 
for the storage changes measured by GRACE or estimated by P — ET. We 354 
further explore this concept, i.e., the negligible amount of fluvial ex- 355 
change between the wetlands and main river channels, in Section 3.3, 356 
below. Instead of river supply, the other potentially significant supply 357 
of water to the wetlands is runoff from the surrounding uplands and 358 
rainfall directly on the wetlands. The P — ET runoff volume anomalies 359 
agree well with the GRACE wetland volume anomalies in terms of tim- 360 
ing and amplitude in region 1 and reasonably well in region 4. In region 361 
2, there is a large discrepancy in their amplitudes in 2003 and 2004. 362 
However, both the GRACE and P — ET anomalies show similar trends 363 
throughout the six years time period. The P — ET runoff volume anom- 364 
alies and the GRACE wetland volume anomalies both reveal a dry 365 


Table 2 t2.1 

Summary of the hydrologic characteristics for the major rivers which have contributing 
areas greater than 10,000 km 2 based on simulation results from the HRR model for the 
period 2003-2008. 


ID 

Study 
region ID 

Contributing area 

(km 2 ) 

Annual discharge Annual discharge 

(m 3 s -1 ) per unit area 

(m 3 s -1 km -2 ) 

- t2.2 
t2.3 

1 

1 

31,270 

429 

0.014 

t2.4 

2 

1 

479,839 

5527 

0.012 

t2.5 

3 

1 

45,642 

675 

0.015 

t2.6 

4 

1 

1,348,434 

18,686 

0.014 

t2.7 

5 

3 

16,073 

209 

0.013 

t2.8 

6 

3 

21,639 

322 

0.015 

t2.9 

8 

3 

39,284 

419 

0.011 

t2.10 

9 

3 

18,809 

343 

0.018 

t2.11 

10 

4 

47,199 

720 

0.015 

t2.12 

11 

4 

453,653 

6587 

0.015 

t2.13 

12 

4 

65,691 

1279 

0.019 

t2.14 

13 

4 

25,752 

551 

0.021 

t2.15 

14 

4 

35,334 

753 

0.021 

t2.16 

15 

4 

136,132 

2187 

0.016 

t2.17 

16 

2 

10,547 

283 

0.027 

t2.18 

17 

2 

14,145 

221 

0.016 

t2.19 

18 

2 

13,713 

190 

0.014 

t2.20 

19 

2 

158,137 

2825 

0.018 

t2.21 





Mean: 0.016 

t2.22 
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Region 1 


Region 2 



Region 3 


CM 


Year 



Region 4 



Year 


Fig. 6. Water volume anomalies of river (blue) and runoff (red). The shading illustrates the range wetland water volume anomalies estimated using CSR, JPL, GFZ, and ITG GRACE 
solutions. The black solid line indicates the mean of the four estimates. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of 
this article.) 


season in 2005 and a rather wet season in late 2007. Essentially, both 
data show somewhat wetter years in 2003 and 2004, dryer years in 
2005 and 2006, and returning to wetter years in 2007 and 2008. In re- 
gion 3, we observe that the timing of volume increases and decreases 
do not generally agree between the wetland and the runoff anomalies, 
although they both show an excessive volume of water in 2007. In sum- 
mary, region 1 annually fills and drains about 20 km 3 to 25 km 3 of water 
each year whereas regions 2, 3 and 4 fill and drain about 10 km 3 to 
20 km 3 . 

The differences observed in regions 2 and 3 may be due to problems 
with satellite rainfall products in the equatorial region. For example, 
Beighley et al. (201 1 ) used three satellite derived precipitation datasets 
(TRMM, CMORPH, PERSIANN) to drive the HRR model throughout the 
Congo Basin. The results, which were compared to historical discharges, 
Envisat altimetry measurements and GRACE water storages, show that 


satellite precipitation products provide unreasonably high rainfall for 381 
specific time periods (e.g., all three in Oct-Nov; only CMORPH and PER- 382 
SIANN in Mar-Apr) in the equatorial regions. These findings are also 383 
consistent with previous studies that found large discrepancies between 384 
gage and satellite precipitation over equatorial regions of Africa (e.g., 385 
McCollum et al., 2000; Nicholson et al., 2003). Although additional re- 386 
search is needed to resolve this issue, one possible cause may be related 387 
to the significant level of lightning activity in the region (Williams and 388 
Satori, 2004). 389 

The rates of wetland filling and draining (Fig. 7) are computed by 390 
taking the temporal derivative of the storage anomalies in Fig. 6 391 
(Alsdorf et al., 2010). Regions 1 through 4 have about ±1000 m 3 s -1 392 
to 2000m 3 s _1 of wetland discharge during flooding and draining. 393 
Summing the maximum and minimum wetland flux rates for all four 394 
regions yields ±6400 m 3 s -1 during flooding and emptying, or ± 16% 395 



Fig. 7. Wetland flow rates (blue, left Y-axis) and P — ET anomalies (green, right Y-axis). (For interpretation of the references to color in this figure legend, the reader is referred to 
the web version of this article.) 
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of the mainstem annually averaged discharge, i.e., 40,000 m 3 s _ 1 at the 
historic Kinshasa gage (Fig. 1 ). Summing the maximum volumes for all 
four regions yields 111 km 3 of water stored and subsequently emptied 
each year from the Congo wetland. This corresponds to only about 8% of 
the total volume of water annually discharged from the Congo. 

The timing of wetland filling (when the flux becomes positive from 
negative) and draining (when the flux becomes negative from positive) 
can be compared with the timing of increasing and decreasing of P — ET 
to examine a temporal connection between them. Note that in Fig. 7 we 
are comparing wetland flow rates derived from GRACE (blue line) to 
changes in P — ET (green line). In regions 1 , 2, and 4, P — ET is always in- 
creasing when the wetland flux rates change from negative to positive, 
i.e., from draining to filling. P — ET therefore always comes before the 
wetland filling, and thus from a temporal perspective, the wetland 
infilling starts with the P — ET runoff from the surrounding uplands. 
On the other hand, when the wetland flux rates switch from positive 
to negative, P — ET is always on the decreasing limb of the annual rain- 
fall. This again is expected where the wetland receives the majority of 
its water from upland runoff. 

3.3. Hydraulic analysis using altimeter measurements in the Congo and 
Amazon basins 

Our interpretation that the Congo wetlands receive the majority of 
their water from upland runoff as opposed to exchange with adjacent 
major tributaries as suggested by comparisons of GRACE anomalies 
with P — ET anomalies, is also supported by Envisat altimeter obser- 
vations. The water elevation changes over the wetland regions, 
which have low topographic relief and higher radar backscatter, are 
generated and compared with the water elevation changes over adja- 
cent river channels (for example, red circles in Fig. 8(a)). The vertical 
datum of both river and wetland water elevations is referenced to the 
Earth Gravitational Model 2008 geoid (EGM08; Pavlis et al., 2008). 
Top panels of Fig. 9(a) and (b) show the surface height profiles 
along the altimeter tracks obtained from several altimeter samplings 
over the red circle regions in Fig. 8(a). We observe fluctuations in the 
water elevations of the Congo mainstem and its adjacent wetlands in 
each altimeter overpass. We then generate water elevation change 
time series by combining successive overpasses. It should be noted 
that the wetland regions closest to the river channels along the tracks 


7 

are selected and compared with the river water fluctuations. As 434 
shown in the bottom panels of Fig. 9(a) and (b), the range in wetland 435 
water levels is small compared to the river. Moreover, the wetland 436 
water topographic elevations are overall between 0.5 and 2.5 m 437 
above the river, i.e., the wetland water levels are always greater 438 
than the river. Clearly, the river cannot flow “uphill” into the wet- 439 
lands. Therefore, the wetlands do not receive water from the adjacent 440 
major tributaries or mainstem Congo River and instead can only sup- 441 
ply water to the rivers. While these observations are necessarily local 442 
to the red-circled areas in Fig. 8(a), they support the previous inter- 443 
pretation which used GRACE observations to suggest that the local 444 
upland runoff is the main source of the Congo wetland water. We 445 
have investigated several more altimeter overpass locations (black- 446 
circled in Fig. 8(a)), where the altimeter footprint allows delineation 447 
of wetlands and rivers, and nearly all locations demonstrate that wet- 448 
land water elevations are consistently higher in elevation than the 449 
adjacent river. It may be argued that the river channel at the location 450 
of the altimeter transect is located further downstream than the wet- 451 
land. In this case, the wetland water levels can always be higher than 452 
the river, and it does not necessarily indicate that the river waters 453 
cannot flow into the wetlands. However, this is true only if there 454 
are abundant floodplain channels that connect the wetlands to the 455 
adjacent river channel. Jung et al. (2010b) highlighted the fundamen- 456 
tal differences in the water level changes between the Amazon and 457 
Congo wetlands, using Interferometric Synthetic Aperture Radar 458 
(InSAR) measurements, due to differences in the connectivity of the 459 
floodplain-river systems. The result suggests that connectivity of the 460 
Congo River to the interfluvial wetland area is limited, compared 461 
with the Amazon. 462 

To further demonstrate that altimeter measurements are a hydraulic 463 
indicator of the direction that water can or cannot flow, we also exam- 464 
ine Envisat altimeter measurements over the Amazon Basin. These 465 
serve as a comparison to the Congo. We generate the time series of 466 
water elevation changes over the Amazon mainstem (Solimoes River) 467 
and its adjacent floodplain (Fig. 9(c)). We note that the water levels at 468 
mid-rising stage in the floodplain are lower than the river, but the low 469 
water levels are almost identical. This implies that the water is flowing 470 
down the hydraulic slope from the river to the floodplain during mid- 471 
rising stage. Moreover, in nearly every year, the river clearly rises before 472 
the floodplain. These timings, in combination with the elevations noted 473 



Fig. 8. Color-coded lines represent Envisat 18-Hz ICE-1 retracked surface heights, referenced to EGM96 geoid, over (a) the Congo Basin from cycle 12 (December 2002), and (b) the 
Amazon Basin from cycle 18 (August 2003). Background is SRTM DEM. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of 
this article.) 
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Fig. 9. Water level fluctuations over the Congo mainstem and its adjacent wetlands along (a) Envisat passes 472 and (b) 887 over the red circle regions in Fig. 8 . Water level variations over 
the Amazon River and its floodplain are shown in (c). They are referenced to the EGM08 geoid. (For interpretation of the references to color in this figure legend, the reader is referred to 
the web version of this article.) 


above, imply that the Amazon river is flowing into the floodplain. We 
also observe that the river has a “shoulder” in the mid-rising times of 
every year where the water level stops rising, or at least slows its rising 
before rapidly rising again. The floodplain shows a similar shoulder. This 
further suggests that the floodplain is responding to the river. These ob- 
servations indicate that the floodplain of the Amazon mainstem derives 
its water more from the river, than from the uplands: a conclusion al- 
ready supported by GRACE observations in Alsdorf et al. (2010). 

4. Conclusions 

The amount of water annually filling and draining the Congo wetland 
is large, but only about one-third of the mainstem Amazon floodplain 
volumes (111 km 3 compared to 285 km 3 ; Alsdorf et al., 2010). Based 
on the amplitude comparison among the water volume changes and 
the timing comparison among their fluxes, we conclude that the local 
upland runoff is the main source of the Congo wetland water, not the 
fluvial process of river-wetland water exchange. Delineating whether 
the water comes from local uplands or from distal places via fluvial 
transport presumably makes a difference in the sediment supplies and 
in the carbon and nutrient exchanges. For example, given the hydraulic 
gradient analysis of Section 3.3, it is unlikely that the Congo wetlands are 
filled with sediments derived from erosion processes in the Congo head- 
waters around the western flanks of the East Africa Rift system. Our 
analysis using altimeter measurements, although they could be local ob- 
servations considering the vast size of the basins, supports our conclu- 
sion, highlighting the difference between the Congo wetland and the 
Amazon floodplain hydraulics. Our finding is in alignment with Jung et 
al. (2010b) which concluded that flow patterns in the Congo are less 
governed by channel connectivity because flooded areas in the Congo 
are broadly distributed and do not have abundant floodplain channels 
as in the Amazon. 

Although we assumed the contribution of soil moisture and ground- 
water variation to the total storage change is negligible compared to 
that of the surface water, further studies are needed to accurately deter- 
mine the portions of the soil moisture and groundwater changes that 
account for the total changes in the water balance, compared to the 
channel and wetland discharges. The HRR hydrologic and hydraulic 
model (Beighley et al, 2009, 2011) and the CaMa-Flood macro-scale 
floodplain model (Yamazaki et al., 2011) can help us determine not 


only those portions, but also simulate the wetland storage changes in 512 
the Congo to compare with our results. 513 

514 

Adler, R. F., Huffman, G. J., Chang, A., Ferraro, R., Xie, P., Janowiak, J., et al. (2003). The 515 
version 2 global precipitation climatology project (GPCP) monthly precipitation 516 
analysis ( 1979-present). Journal of Hydrometeorology, 4, 1147-1167. 517 

Alsdorf, D., Han, S. -C., Bates, P., & Melack, J. (2010). Seasonal water storage on the Amazon 518 
floodplain measured from satellites. Remote Sensing of Environment, 1 14, 2448-2456. 519 
Beighley, R. E., Eggert, K. G., Dunne, T., He, Y., Gummadi, V., & Verdin, K. L. (2009). Sim- 520 
ulating hydrologic and hydraulic processes throughout the Amazon River Basin. 521 
Hydrological Processes, 23, 1221-1235. 522 

Beighley, R. E., Ray, R. L., He, Y., Lee, H., Schaller, L., Durand, M., et al. (201 1 ). Comparing 523 
satellite derived precipitation datasets using the Hillslope River Routing (HRR) 524 
model in the Congo River Basin. Hydrological Processes, doi:10. 1002/hyp. 8045 525 
(published online). 526 

Bettadpur, S. (2007). CSR Level-2 processing standards document for product Release 527 
04. GRACE 327-742, Center for Space Research, University of Texas at Austin. 528 

Bwangoy, J. -R. B., Hansen, M. C., Roy, D. P., de Grandi, G., & Justice, C. O. (2009). Wet- 529 
land mapping in the Congo Basin using optical and radar remotely sensed data and 530 
derived topographical indices. Remote Sensing of Environment, 1 14, 73-86. 531 

Duan, X. J., Guo, J. Y., Shum, C., 8 * van der Wal, W. (2009). On the postprocessing removal of 532 
correlated errors in GRACE temporal gravity field solutions. Journal of Geodesy, 83, 533 
1095-1106. 534 

Guo, J. Y., Duan, X. J., & Shum, C. (2010). Non-isotropic filtering and leakage reduction 535 
for determining mass changes over land and ocean using GRACE data. Geophysical 536 
Journal International, doi:10.1111/j.l365- 246X.201 0.04534.x. 537 

Han, S. -C., Shum, C., Jekeli, C., & Alsdorf, D. (2005). Improved estimation of con- 538 
tinental water storage change from GRACE. Geophysical Research Letters, 32, 539 
doi:10.1029/2005GL022382. 540 

Hansen, M. C., Roy, D. P., Lindquis, E., Adusel, B., Justice, C. O., & Altstatt, A. (2008). A method 541 
for integrating MODIS and Landsat data for systematic monitoring of forest cover and 542 
change in the Congo Basin. Remote Sensing of Environment, 1 12, 2495-2513. 543 

Jung, H., Alsdorf, D., Lee, H., Trigg, M., 82 Fewtrell, T. (2010). Hydrogeomorphic flood 544 
classification and hydrodynamic modeling of the Congo interfluvial wetlands. 545 
AGU fall meeting abstract. 546 

Jung, H., Hamski, J., Durand, M., Alsdorf, D., Hossain, F., Lee, H., et al. (2010). Character- 547 
ization of complex fluvial systems using remote sensing of spatial and temporal 548 
water level variations in the Amazon, Congo, and Brahmaputra Rivers. Earth Sur- 549 
face Processes and Landforms, 35, 294-304. 550 

Laporte, N., Justice, C., & Kendall, J. (1995). Mapping the dense humid forest of Cameroon 551 
and Zaire using AVHRR satellite data. International Journal of Remote Sensing, 16, 552 
1127-1145. 553 

Laraque, A., Mahe, G., Orange, D., & Marieu, B. (2001 ). Spatiotemporal variations in hy- 554 
drological regimes within Central Africa during the XXth century .Journal ofHydrol- 555 
ogy, 245, 104-117. 556 

Lee, H., Durand, M., Jung, H., Alsdorf, D., Shum, C., & Sheng, Y. (2010). Characterization 557 
of surface water storage changes in Arctic lakes using simulated SWOT measure- 558 
ments. International Journal of Remote Sensing, 31, 3931-3953. 559 

Lehner, B., Verdin, K., & Jarvis, A. (2008). New global hydrography derived from space- 560 
borne elevation data. Eos, 89, 93-94. 561 


References 


Please cite this article as: Lee, H., et al., Characterization of terrestrial water dynamics in the Congo Basin using GRACE and satellite radar altimetry, 
Remote Sensing of Environment (201 1 ), doi: 10.1 01 6/j.rse.201 1.08.01 5 






562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

589 


ARTICLE IN PRESS 


H. Lee et al. / Remote Sensing of Environment xxx (201 1 ) xxx-xxx 9 


McCollum, J. R., Gruber, A., & Ba, M. B. (2000). Discrepancy between gauge and satellite es- 
timates of rainfall in equatorial Africa .Journal of Applied Meteorology, 39, 666-679. 
Melack, J. M., & Engle, D. (2009). An organic carbon budget for an Amazon floodplain 
lake. Verhandlungen — Internationale Vereinigung Limnologie, 30, 1179-1182. 
Nicholson, S. E., Some, B., McCollum, J., Nelkin, E., Klotter, D., Berte, Y., et al. (2003). Val- 
idation of TRMM and other rainfall estimates with a high-density gauge dataset for 
West Africa. Part I: Validation of GPCP rainfall product and pre-TRMM satellite and 
blended products. Journal of Applied Meteorology, 42, 1337-1354. 

Pavlis, N. K., Holmes, S. A., Kenyon, S. C., & Factor, J. K. (2008, April). An earth gravita- 
tional model to degree 2160: EGM2008. Vienna, Austria: EGU-2008 (also see). 
http://earth-info.nga.mil/GrandG/wgs84/gravitymod/egm2008/index.html 
Swenson, S., & Wahr, J. (2006). Post-processing removal of correlated errors in GRACE 
data. Geophysical Research Letters, 33, L08402, doi:10.1029/2005GL025285. 


Tapley, B. D., Bettadpur, S., Ries, J., Thompson, P., & Watkins, M. (2004). GRACE mea- 575 
surements of mass variability in the Earth system. Science, 305, 503-505. 576 

Verdin, K. L., & Verdin, J. P. ( 1 999 ) . A topological system for delineation and codification 5 77 
of the earth's river basins .Journal of Hydrology, 218, 1-12. 578 

Wahr, J., Molenaar, M., & Bryan, F. (1998). Time variablility of the earth's gravity field: 579 
Hydrological and oceanic effects and their possible detection using GRACE .Journal 580 
of Geophysical Research, 103, 30205-30229. 581 

Williams, E. R., & Satori, G. (2004). Lightning, thermodynamics and hydrological com- 582 
parison of the two tropical continental chimneys. Journal of Atmospheric and Solar- 583 
Terrestrial Physics, 66, 1213-1231, doi: 10.101 6/j.jastp. 2004.05.015. 584 

Yamazaki, D., Kanae, S., Kim, H., & Old, T. (2011). A physically based description of 585 
floodplain inundation dynamics in a global river routing model. Water Resources 586 
Research, 47, W04501, doi:10.1029/2010WR009726. 587 

588 


Please cite this article as: Lee, H., et al., Characterization of terrestrial water dynamics in the Congo Basin using GRACE and satellite radar altimetry, 
Remote Sensing of Environment (201 1 ), doi:10.1016/j.rse.201 1.08.015 




